Are non-equilibrium Bose-Einstein condensates superfiuid? 
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We theoretically study the superfluidity properties of a non-equilibrium Bose-Einstein condensate 
of exciton-polaritons in a semiconductor microcavity. The dynamics of the condensate is described 
at mean-field level in terms of a modified Complex Ginzburg Landau equation. A generalized 
Landau criterion is formulated which estimates the onset of the drag force on a small moving defect. 
Metastability of supercurrents in multiply connected geometries persists up to higher flow speeds. 
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Superfluidity is among the most remarkable conse- 
quences of macroscopic quantum coherence in condensed 
matter systems and manifests itself in a number of fasci- 
nating effects [3, 0] ■ A unified description of these phe- 
nomena is obtained in the framework of the so-called two- 
fluid hydrodynamics, in which the macroscopic conden- 
sate wavefunction adds up to the standard hydrodynamic 
variables [3[ . The phenomenon of macroscopic coherence 
is not restricted to systems at (or close to) thermodynam- 
ical equilibrium such as liquid Helium, ultracold atomic 
gases, or superconducting materials, but has been ob- 
served also in systems far from thermodynamical equi- 
librium, whose state is determined by a dynamical bal- 
ance of driving and losses. Most remarkable examples 
are lasers and, more recently, Bose-Einstein condensates 
of magnons in magnetic solids Hand exciton-polaritons 
in semiconductor microcavities [5|, |fj( . In particular, the 
issue of superfluidity in this latter system has attracted 
a significant interest from both the theoretical @, H[ and 
experimental [9Mll| points of view. 

Recent experiments with resonantly pumped polari- 
ton condensates [t| have demonstrated superfluidity as a 
dramatic reduction in the intensity of resonant Rayleigh 
scattering as originally predicted in ML The situation is 
less clear in the case of non-resonant [5( or parametrical 
(OPO) [Io| pumping schemes: recent experiments (l(i| 
have observed propagation of polariton bullets without 
apparent friction, which is in contrast with the predic- 
tions of a naive Landau criterion based on the elementary 
excitation spectrum predicted in [l2l ]. Another aspect 
of superfluidity, namely metastability of supercurrent in 
multiply-connected geometries was investigated theoret- 
ically in Q and experimentally confirmed in [ill ]. The 
present paper reports a comprehensive theoretical inves- 
tigation of the meaning of superfluidity for polariton con- 
densates under a non-resonant pumping. Emphasis will 
be given to the novel features that originate from their 
non-equilibrium character. 

At the mean-field level, the condensate dynamics can 
be described in terms of the so-called Gross-Pitaevskii 
equation (GPE) 0], which was recently generalized to 
non-equilibrium condensates by including the effect of 



pumping and losses [13|, [14[. This mean- field descrip- 



tion has been able to explain a number of experimen- 
tal observations on polariton condensates, e.g. the ring- 
shaped momentum distribution of spatially narrow con- 
densates la . Ilal . the synchronization/desynchronization 
transition [la ], the spontaneous appearance of vor- 
tices Nonetheless, the implicit assumption that the 
pumping mechanism is not frequency-selective can lead 
to unphysical predictions, e.g. that in a spatially ho- 
mogeneous or ring-like geometry condensation is equally 
likely to occur in any momentum state. Kinetic calcu- 
lations [181 ] have pointed out the significant energy de- 
pendence of the polariton-polariton scattering processes 
that are responsible for replenishing the condensate. In- 
cluding this feature as an energy-dependent amplification 
mechanism turned out to be crucial in order to extract 
physically meaningful predictions for the condensate fluc- 
tuations in the Wigner Monte Carlo simulations of [lil ]. 

A simplest generalization of the GPE to include 
frequency-dependent pumping has the form: 
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The energy zero has been set for convenience at the bot- 
tom of the lower polariton branch; the efficiency of am- 
plification (proportional to the pumping strength P) de- 
creases to zero a frequency interval Qk above it. As- 
suming a linear form of the frequency dependence of the 
amplification, the generalized GPE (JTJ) maintains a tem- 
porally local form. The others terms describing gain sat- 
uration (r), losses (7), polariton mass (m), polariton- 
polariton interactions (<?), external potential (V e xt) have 
the same meaning as in previous works [HI, [3] ■ 

We first consider the evolution of the system starting 
from an initial state with no condensate ip = 0. If the 
strength of pumping is enough to overcome losses P > 7, 
the -0 = state is dynamically unstable against the cre- 
ation of a finite condensate amplitude in all the low- 
momentum modes for which hk 2 /2m < fix (1 — l/P)- 
The rate of this instability is maximum at k = and 
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decreases for increasing k. This fact is in agreement with 
the experimental observation of condensation naturally 
occurring around k = as soon as the sample is suffi- 
ciently large and free from disorder 

In spite of this natural preference, condensation can be 
forced to occur in finite momentum state by seeding the 
system with a short resonant pulse at the desired k c as 
proposed and numerically assessed in Ref. A related 
configuration was experimentally demonstrated for the 
OPO pump scheme in The effect of the frequency- 
dependent pumping is then visible in the state equation 
relating the pumping strength P to the condensate den- 
sity n c = \ip\ 2 in the stationary state: 
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As a consequence of the frequency dependent pumping, 
the interaction-induced blue-shift of the condensate fre- 
quency oj c — fc 2 /2m + g \ip | 2 is responsible for a slower 
increase of density with pumping strength P. 
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FIG. 1: Elementary excitation spectrum of a moving, spa- 
tially homogeneous, non-equilibrium condensate as described 
by the linearized CGLE around the steady-state solution. 
Parameters: fe c /fc 7 = 4, r/f = 1, g/j = 1, Ok/7 = 50, 
n cdll — 0-1 ( c >d) and n c g/"f = 1 (e,f). Inset: magnified 
view of (d). Panels (a,b): same as panels (e,f) in the limit 
of a frequency-independent pumping Qk = oo. Momenta are 
measured in units of k~. 



The dynamical stability of a moving condensate is to 
be assessed by linearization of the CGLE around the 
stationary solution. Stability requires that the imagi- 
nary part of the frequency is negative for all Bogoliubov 
modes. Examples of the Bogoliubov dispersion uJBogil) 
for moving condensates are shown in Fig[TJ In the left- 
most panels (a,b), the limit fix — > oo of a negligible 
frequency-dependence of pumping is considered. Apart 
for the global Doppler tilting of the real part due to the 
finite condensate velocity, the plotted dispersion fully re- 
covers the diffusive character first discussed in 

Hi-El]: a 



flat dispersion around the condensate wavevector and an 
imaginary part quadratically growing in q. In this limit, 
stable condensates exist for any value of k c . 

The effect of a frequency-selective pumping is ad- 
dressed in the other panels. The central panels (c,d) refer 
to the case of a pumping intensity just above the thresh- 
old. In this regime, the condensate density is very small 
and the real part of the Bogoliubov mode frequency re- 
duces to the single-particle one. On the other hand, the 
characteristic damping rate of density fluctuations (the 
gapped mode at q — 0) and of the high-momentum modes 
is suppressed by a critical slowing down phenomenon [p| . 
As a result, the low energy modes around q ~ — k c be- 
come dynamically unstable. Eventually, this mechanism 
leads to the disappearence of the original condensate and 
the formation of another -stable- condensate in a lower 
momentum state. 




FIG. 2: Growth rate (in units of 7) of maximally unstable 
mode as a function of condensate momentum k c and density 
n c . The dashed line indicates the boundary of the stability 
region. Same system parameters as in Figfjjc-f). 

The right panels (e,f) refer to the case of a stronger 
pumping well above the threshold. In this case, the 
damping rate of all modes other than the Goldstone one is 
comparable to the polariton lifetime 7. In this case, the 
frequency-dependence of pumping is not able to desta- 
bilize the moving condensate. The stability domain as 
a function of the momentum k c and the density n c of 
the condensate is summarized in Fig{2j as the density 
is increased, stable condensates survive up to larger mo- 
menta. 

It is of crucial importance to note that this instability 
mechanism has no direct counterpart in standard equi- 
librium superfluids and is effective even in the absence 
of any defect potential. Of course, the present theory is 
directly applicable only to infinite, homogeneous conden- 
sates or in multiply connected ring-shaped geometries. In 
finite systems, assessing the stability of moving conden- 
sates would be made harder by the flow of condensate 
polaritons outside the pump regions (ljj ]. 

In the presence of weak defects, the frictionless flow of 
equilibrium condensates is limited to flow speeds below 
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the Landau critical velocity v c — min^ {Re[a;^ og (A:)]/A;}, 
where oj° Bog (k) is the real part of the Bogoliubov dis- 
persion for a condensate at rest. For dilute systems, v c 
coincides with the sound velocity c s — -J gn c /m [l|. As 
a consequence of the diffusive character of the Goldstone 
mode, a naive application of this definition to the Bo- 
goliubov dis per sion of non-equilibrium condensates an- 
ticipated in [12h14| | gives a vanishing prediction for the 
critical velocity, v c = 0: Bogoliubov waves are emitted in 
a moving condensate hitting a (weak) defect at any value 
of the condensate speed. 
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FIG. 3: Density perturbation created in a moving condensate 
by a stationary weak defect for three different values of the 
condensate velocity v/c a — 1.5, 1, 0.4 across the critical value 
for superfluidity. Parameters: n c g/^f = n c r/7 = 1, Qfc/'y = 
50. 



A complete picture of non-equilibrium condensates re- 
quires taking into account the non-trivial dynamics of 
the imaginary part of the Bogoliubov dispersion. Nu- 
merical plots of the density perturbation induced by a 
single weak stationary defect in a (dynamically stable) 
moving condensate can be calculated from the general- 
ized GPE (TTJ in the presence of a suitable defect potential 
V ex t(r) and are shown in Fig J3] for different values of the 
condensate velocity. In contrast to the predictions of the 
naive Landau criterion stated above, the induced per- 
turbation closely resembles the one induced in an equi- 
librium condensate described by the standard GPE [221 ]: 
at high speeds [panel (a)], the defect creates a series of 
parabolic fringes that propagate away from the defect; 
at low speeds [panel (c)], the propagating fringes are re- 
placed by a localized perturbation in the vicinity of the 
defect. Remarkably, the characteristic speed at which 
the fringes disappear is of the order of the equilibrium 
sound speed c s = ^gn c /m, but does not correspond to 
any noticeable feature in the real part of the Bogoliubov 
dispersion of Fig|TJa,c,e). 

As it was pointed out in a different context in [20| , the 
emission pattern by a localized monochromatic source is 
better understood in terms of the (complex) wavevector 
of the emitted wave at the given (real) frequency. As one 
can see in Fig. 0] (b) , the wavevector k of the Bogoliubov 
wave emitted at zero frequency by the static defect starts 
having a finite real part only after a branching point: 
extended oscillations in the density are observed as soon 
as the real part of k exceeds the imaginary part. 
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FIG. 4: Panel (a): Force exerted on a weak stationary defect 
by a moving condensate as a function of the condensate veloc- 
ity for different values of rn c /"/ = 0.1, 1,2 (blue solid, green 
dashed, red dotted lines), in units of F 7 = ft 2 fc 7 /m. For com- 
parison, the same force on an equilibrium condensate is shown 
as a thin solid black line. Panel (b,c): Real and imaginary 
parts of the complex wavevector k of the Bogoliubov emission 
in the upstream direction as a function of the condensate ve- 
locity, for rn c /7 = gn c /'y — 1. Fhe other parameters are the 
same as in Fig. [3] 



This behavior is reflected in the drag force F exerted 
by the moving condensate onto the defect as a function 
of the condensate velocity [FigH^a)]. In terms of the 
perturbed density profile n(r\ the drag force is given by 
F = — J d 2 r n(r) V r Vd e /(r) [2lJ. As a consequence of the 
finite lifetime of the Bogoliubov modes, the drag force 
has a non-vanishing value at all v. In addition, the drag 
force shows a marked threshold at a velocity value that 
closely corresponds to the onset of fringes in the density 
profile: the lower the value of the n c r/j parameter, the 
sharper the threshold. On the other hand, the role of fix 
on this effect is minor. Even though no experiment has 
so far investigated the superfluid properties of polariton 
condensates under non-resonant pumping, it is likely that 
this generalized form of the Landau criterion provides at 
least a partial explaination of the recent ex peri mental 
observation of superfluidity in a OPO regime [id ]. 

Supersonic flows in equilibrium condensates in ring- 
shaped geometries are generally strongly sensitive to the 
presence of defects: as soon as nodes appear in the con- 
densate wavefunction, the topological stability of the su- 
percurrent state is broken, which leads to a rapid slow 
down of the condensate motion. The finite damping 
rate of excitations in polariton condensates introduces 
a substantial modification to this picture: the perturba- 
tion created by each defect is not able to propagate on 
long distances, but rather remains localized in space on 
a length scale inversely proportional to Im[fc]. As a re- 
sult, we expect it is much harder to break the topological 
stability of the supercurrents. 

This physical picture is confirmed by numerical simu- 
lations of the time evolution under the generalized GPE 
([I]) starting from an initial condition with a supersonic 
flow. Periodic boundary conditions are assumed. For a 
low defect density, the supercurrent state is maintained 
for very long times with no sign of decay. The char- 
acteristic Cerenkov-like density patterns in the vicinity 



4 




1 



-2 2 

k /k 





a 



2 

k /k 

x y 



modes and preventing the momentum distribution from 
spreading all over the fc-space. 

In conclusion, we have investigated how the non- 
equilibrium nature of exciton-polariton condensates af- 
fects their superfluidity properties. A mean-field model 
based on a generalized Gross-Pitaevskii equation includ- 
ing a frequency-dependent pumping is developed. Dif- 
ferent aspects of superfluidity have been considered. A 
non-equilibrium version of the Landau critical speed for 
the onset of drag force is formulated and metastability 
of supercurrents is shown to persist even at speeds well 
above the critical speed. 

We are indebted to Vincenzo Savona, Cristiano Ciuti 
and Davide Sarchi for continuous enlightening exchanges. 
Stimulating discussions with J. Keeling, C. Menotti, F. 
Piazza, D. Sanvitto, A. Smerzi are acknowledged. 



FIG. 5: Real (linear color scale) and momentum space (loga- 
rithmic color scale) densities of a non-equilibrium condensate 
after a temporal evolution of "/t = 300 for an initial momen- 
tum k/kj = 2 and two different densities of defects. The dots 
in the real space panels indicate the position of the defects. 
The square in the momentum space panels indicates the ini- 
tial momentum of the condensate. Parameters are the same 
as in Fig. 2] except for Q.k/i = 10. 



of each defect are spatially separated and do not inter- 
fere [FigJS^a)]. The momentum distribution in /c-space 
[FigJSJb)] shows the condensate peak right at the initial 
momentum state and a much fainter resonant Rayleigh 
scattering ring [7[ . The situation is different for a higher 
density of defects. In this case, the real space density per- 
turbations created by the different defects substantially 
overlap with each other [FigEJc)]. As a result, interfer- 
ence effects are more likely to create nodes in the con- 
densate wavefunction and therefore to trigger dissipation 
of the supercurrent. The occurrence of such a process is 
apparent in FigEJd) as a much reduced late-time value 
of the condensate momentum. 

This fact illustrates another important difference with 
respect to standard, equilibrium condensates: in that 
case, the dissipation of a supercurrent leads to a signifi- 
cant heating and reduction of the condensate fraction. In 
the present case, the condensate turns out to be almost 
rigidly transferred to a lower momentum state; the (co- 
herent) momentum broadening that is visible in FigJSfd) 
is an effect of the defects. Remarkably, the frequency de- 
pendence of pumping in the generalized CGLE (p} has a 
crucial role in concentrating the population in low-energy 
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